Phase diagrams of the metallic zigzag carbon nanotube 
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We investigate a metallic zigzag carbon nanotube by means of a Hubbard model which includes 
both on-site and nearest neighbour interactions. Assuming weak interactions, a renormalization 
group analysis of the equivalent two-leg ladder followed by bosonization and refermionization results 
in a Gross-Neveu model with an enlarged symmetry relative to the original Hamiltonian. For the 
undoped case the symmetry of the Gross-Neveu model is SO (8), but for the doped case the particle- 
hole symmetry is broken and the symmetry reduces to SO(6). Four ground state phases are found 
in the undoped carbon nanotube with repulsive interactions, a d-wave Mott insulator, an s-wave 
Mott insulator, a p-density wave and a charge density wave. The doped case has two ground state 
phases, a d-wave superconductor and a phase where a p-density wave and a charge density wave 
co-exist. We also explore the global phase diagram with a general interaction profile and find several 
additional states, including a chiral current phase where current flows around the nanotube along 
the zigzag bonds. 

PACS numbers: 61.50.Ah, 73.22.-f, 71.10.Fd, 71.10.Hf 
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I. INTRODUCTION 

Carbon nanotubes (CNT) are long, thin tubes con- 
structed from sheets of graphite. Because of their 
many novel properties CNT have numerous poten- 
tial applications in material science optics^ 5 - and 
electronics jSiLSi while also contributing to our knowledge 
of fundamental physics i 9 ' 10 ' 11 They are extremely strong, 
owing to being constructed from sp 2 bonds, and yet their 
low density makes them extremely light. CNT have the 
highest tensile strength and elastic modulus of any known 
materia l 12 ' 13 but this impressive strength is only applica- 
ble to forces which stretch the nanotube as their hollow 
structure means they readily become distorted under tor- 
sion, compression or bending^ The electrical properties 
of a CNT depend on its structure. Single walled CNT are 
generally classified as one of three types: zigzag, arm- 
chair or chiral, and these classifications are defined by 
the orientation of the graphene lattice about the tube. 
From band structure calculations, armchair CNT are al- 
ways metallic with current densities which may possibly 
exceed silver and copper ji 5 - while zigzag and chiral CNT 
can be metallic, semiconducting or insulating depending 
on the width and helicity of the tube 

The simple band structure calculations used to deter- 
mine whether or not a CNT is metallic are performed 
in the weakly interacting limit. In general this limit is 
not applicable to CNT as they tend to have long-range 
Coulomb interactions which are not small, though one 
can justify considering only weak short-range interactions 
if the CNT is screened. Screening can be achieved either 
by arranging several nanotubes in an array or ropc-like 
structure ; 21 ' 22 or by placing a single nanotube close to 
a conducting platej 23 ' 24 Studies on armchair CNT with 
long-range interactions have shown that the doped CNT 
ground state is a metallic Tomanaga-Luttinger liqui d 25 ' 26 
while the ground state of the undoped CNT has a number 



of possible phases including Mott insulators and density 
Some of these phases have also been found 



waves 



27,28 



in doped and undoped armchair CNT with short-range 



interactions 



29,30 



Here we consider a metallic zigzag CNT which we 
assume to be screened so that only on-site and near- 
est neighbour interactions need to be considered. Our 
main purpose is to determine which ground state phases 
can be supported by a metallic zigzag CNT with short- 
range interactions, and to describe the phase transitions 
between these phases. The Hubbard Hamiltonian of a 
metallic CNT may be mapped onto the well known Hub- 
bard Hamiltonian of a two-leg ladder ; 31 ' 32 ' 33 though the 
nature of the two-leg ladder is dependent on the chirality 
of the nanotube. The armchair CNT maps onto a fairly 
standard form of the two-leg ladder where hopping be- 
tween any two adjacent lattice sites along either leg is al- 
ways the same, as is hopping along any rung i.e., between 
legs. So, it is not hard to determine the behaviour of 
an armchair CNT directly from known results of two-leg 
laddersj 29 ' 30 ^ 3 ^ 3 ^^ 7 . In contrast, the metallic zigzag 
CNT maps onto an unusual type of two-leg ladder in 
which the hopping part describes two chains with hop- 
ping strength alternating between lattice sites, but no 
hopping between legs. The two chains however cannot 
be described as independent as they influence each other 
through Coulomb interactions. 

The nearest neighbour Coulomb interactions in an 
armchair CNT map rather simply onto a two-leg lad- 
der, acting between nearest neighbours either along the 
rungs or along the legsi^ 3 . The situation is quite differ- 
ent in the two-leg ladder equivalent of a zigzag CNT. In 
this case the interactions in the ladder act either between 
nearest neighbours along the legs or between next near- 
est neighbours on different legs. Initially we hoped that 
these more complex interactions would have an interest- 
ing effect on the phase diagram, possibly allowing some 
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unusual phases with broken time-reverasl symmetry. For 
instance, it is interesting to explore whether a staggered 
flux phase exists under physically possible conditions. As 
we shall show later, this is not the case and the physi- 
cally possible phases of a zigzag CNT are similar to those 
found in an armchair CNT. 

We determine that our metallic CNT can support at 
least seven different phases when there is no doping, 
though some of these phases are only found when at- 
tractive interactions are permitted. Four phases are Mott 
insulators, two with ci-wave symmetry which we name D- 
Mott and D'-Mott, and two with s-wave symmetry which 
we name S-Mott and S'-Mott. These four Mott phases 
are analogous to those found in a standard two-leg ladder 
and so we have adopted the same naming convention.— 
The oder parameters of currents and/or bond hopping 
are always zero in the Mott insulator phases, and the 
average electron density is always one electron per site. 
These phases are classified as d-wave or s-wave based on 
the nature of the pairing correlations. Illustrations of all 
four Mott states in both the honeycomb lattice of the 
zigzag CNT and the equivalent two-leg ladder are shown 
in Figs. QJa),(b),(c) and (d). In these figures the circles 
represent an s-wave pairing of two electrons with oppo- 
site spin. One pair must occupy one of two possible sites 
but which site is chosen is completely random, ensuring 
an average electron density of one electron per site. In 
the S-Mott the two possible positions for one electron 
pair are nearest neighbours so we draw circles between 
all nearest neighbours. In the S'-Mott the two possible 
positions for one pair are next-nearest neighbours so the 
circles are drawn between next- nearest neighbours. The 
d-wave pairing is represented by ellipses. In this case 
paired electron occupy different sites, though the posi- 
tions of the different spins are not ordered but random. 
The paired electrons are nearest neighbours in the D- 
Mott so the ellipses are drawn along the bonds joining 
nearest neighbours. In the D'-Mott the pairing is be- 
tween next-nearest neighbours so the ellipses (distorted 
for clarity in the two-leg ladder case) join next-nearest 
neighbour sites. 

The remaining three phases are density waves which 
all exhibit a broken Z2 symmetry. The charge density 
wave (CDW) has broken particle-hole symmetry so there 
is a regular variation in the electron distribution along 
the lattice as shown in Fig. [He). As before the cir- 
cles represent a pair of opposite spin electrons. A major 
difference between the CDW and the two s-wave Mott 
states is that positions of the pairs in the CDW is not 
random. The p-density wave (PDW) is equivalent to a 
spin-Peierls state where dimers form between neighbour- 
ing sites as shown in Fig. QJf). The thickness of the 
dashed and solid black lines indicates the magnitude of 
the kinetic energy exchange between sites, though these 
two types lines have opposite signs. The chiral current 
phase (CCP) describes a state in which a current cir- 
culates around the nanotube, flowing between nearest 
neighbours. This is represented in Fig. []Jg) with the 
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FIG. 1: Eight phases of the zigzag CNT and the equivalent 
two- leg ladder in the strong coupling limit. The symbols are 
defined in the text. 



arrows describing the current. In the equivalent two-leg 
ladder all currents vanish because the two-leg ladder lies 
along the longitudinal axis of CNT and the net current 
along the longitudinal axis of the CNT is zero. There 
is fourth density wave phase which is theoretically possi- 
bly, though we do not find it in our phase diagram. This 
phase is a /-density wave (FDW) which has a circulating 
current which flows between next-nearest neighbours, as 
shown in Fig. [ljh) . The current in the equivalent two-leg 
ladder vanishes, again because the net current along the 
CNT longitudinal axis is zero. 

The eight phases are essentially equivalent to phases 
found in the armchair CNT (though the phase diagram 
is not identical), and comparable to phases found in the 
two-leg ladder, with the exception of the CCP. When an 
armchair CNT is in a CCP the equivalent two-leg ladder 
is in a staggered flux (SF) phase in which current flows 
around plaquettes between nearest neighbours, with the 
direction of the current being opposite in neighbouring 
plaquettes.— The two-leg ladder equivalent of the zigzag 
CNT's CCP resembles a SF phase, yet it cannot be a true 
SF phase as this ladder does not have standard plaquettes 
around which a non-zero current may flow. 

If a two-leg ladder or CNT is lightly doped away 
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from half-filling different phases emerge, yet they are 
still closely related to the phases found in the half-filled 
case. Two phases are superconducting, one being d-wave 
(D-SC) and the other s-wave (S-SC). The D-SC can be 
thought of as a merging of the D-Mott and the D'-Mott 
insulator phases found in the undoped case. Similarly, 
the S-SC is a merging of the undoped S-Mott and S'- 
Mott phases. The doped two-leg ladder and CNT have 
two density wave phases. One is a combination of the 
undoped CDW and PDW phases. For simplicity we will 
refer to this doped phase as a CDW. The other density 
wave is a combination of the CCP and the FDW in the 
CNT (or a SF and FDW in the two- leg ladder), and we 
shall refer to this phase as a CCP. The general appear- 
ance of the phase diagrams of the doped and undoped 
cases are quite similar, although, in the undoped cases 
the Mott insulator phases tend to dominate most of the 
phase diagram, while on doping the density wave regions 
are significantly enhanced. 

One fascinating characteristic of undoped two-leg lad- 
ders, also shown to exist in CNT, is dynamical symmetry 
enlargement (DSE). In the two-leg ladder the symme- 
try of the original Hubbard Hamiltonian is U(l)xSU(2) 
but the effective Hamiltonian obtained after solving the 
renormalization group (RG) equations is equivalent to a 
Gross-Neveu (GN) model with a SO(8) symmetry, which 
is far larger than the original symmetry. 29 All phases 
exhibit DSE, yet they do not all share the same SO (8) 
symmetry. A doped two-leg ladder can still be mapped 
onto a GN model with DSE, but with the smaller SO (6) 
symmetry— It is important to note that DSE is a con- 
sequence of weak coupling and begins to break down as 
the coupling increases^ Furthermore, the SO (8) symme- 
try is not particularly realistic as it depends on an equal 
charge and spin gap, while in reality the charge gap is 
substantially larger than the spin gap. In contrast, the 
SO (6) symmetry is fairly realistic as doping removes the 
charge gap but has no effect on the spin gap. 

In Sec. |H] we show how a Hubbard model describing 
a metallic zigzag CNT at half-filling with both on-site 
and nearest neighbour interactions may be mapped onto 
a two-leg ladder. The Hamiltonian is then bosonized. 
In Sec. IIIII we obtain all phases and phase transitions 
by solving the RG equations under different initial con- 
ditions and substituting the results into the bosonized 
Hamiltonian. We then show how the effective Hamilto- 
nian can be represented as a SO (8) GN model with en- 
larged symmetry. We use a variety of order parameters 
such as the current and the electron density to determine 
the characteristics of each phase. We discuss the nature 
of the phase transitions and show that they can mostly 
be classified as either Gaussian or Ising. In Sec. IIVI 
we consider the case of a doped CNT. The Hamiltonian 
is similar to the undoped Hamiltonian, except for the 
absence of Umklapp interactions. Using a RG analysis 
and bosonization we find the doped nanotube's ground 
state phases and the phase transitions between them and 
we show that all phases map to an effective SO(6) GN 



Hamiltonian. Finally, in Sec. |V]we discuss some general 
properties of the SO(2N) GN model. 

II. THE MODEL 

In this section we show how the zigzag CNT Hubbard 
model with on-site and nearest neighbour interactions 
may be mapped onto a two-leg ladder Hubbard model. 
After some standard approximations we bosonize the 
Hamiltonian. These derivations are not new , 29 ' 33 except 
for the addition of the nearest neighbour interactions, 
and therefore our derivation will be rather brief. 

A carbon graphene lattice may be represented by two 
regular triangular sublattices, offset by d = a(0, — I/V3) 
and with sublattice basis vectors a± = a(±l/2, v / 3/2), 
where a is the sublattice lattice constant, as illustrated 
in Fig. [2 The Hubbard model hopping Hamiltonian of 
such a carbon lattice is^ 3 - 

Ho = -t £ [cL(r)c 2 a(r + d + )+c t lQ (r)c 2Q (r + d_)] 

r— R..Q 

"*-L £ [c\jr)c 2a (r + d) + h.c] (1) 

r— R..Q 

where d± = &± + d, R = n + & + + ?i_a_ with integral 
n± describes a lattice vector in the first sublattice, t± is 
the hopping strength along the vertical (i.e., y direction) 
bond and t is the hopping strength along the other two 
bonds. In all calculations we define t = t±. The annihi- 
lation operator Ci a describes the destruction of a fermion 
with spin a in the ith sublattice. The on-site interaction 
Hamiltonian is 

Hu = U ^2 : n lT n a : (2) 

where rii a — c\ a Ci a , every lattice site is given by R' = 
?i+ a + + n_a_ + n^d with integral n±^, and U is the on- 
site interaction strength. Similarly, the nearest neighbour 
interaction Hamiltonian is 

Hy =V £ [n la (r)n 2 p{r + d+) + n la (r)n 2 p(r + d_)] 

oftrER 

+ V± £ ni Q (r)n 2/3 (r + d) (3) 

where V± is the nearest neighbour interaction strength 
across the vertical bond and V is the nearest neighbour 
interaction strength across the other two bonds. We will 
only consider \V\, |Vj_| < \U\. 

A single walled CNT is formed by making a cylinder 
out of any graphene lattice, such as the one shown in 
Fig. [2] If this lattice is rolled along a horizontal axis so 
that the top is joined to the bottom of the lattice (while 
retaining the regular hexagonal structure) we obtain an 
armchair CNT. If instead we roll this lattice along a ver- 
tical axis, joining the left and right sides, we obtain a 
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FIG. 2: The graphene lattice with the two triangular sublat- 
tices shown as red squares and green triangles. 



zigzag CNT. A chiral CNT is any other cylinder which 
can be created from a graphene lattice which is neither 
armchair nor zigzag. Here we will only consider a metal- 
lic zigzag CNT. To determine when the zigzag CNT is 
metallic we derive the band structure in the weak cou- 
pling limit. 

The y-axis of the zigzag CNT is defined to be in the 
longitudinal direction of the nanotube and the x-axis is 
around the nanotube in the transverse plane. As the 
momentum must be quantized in the x direction 



k x — 



2itp 



p = 0,±l,...,±(N x /2), 



(4) 



where N x is the number of lattice sites around the cir- 
cumference of the CNT and the circumference is aN x . 
The CNT can be defined as metallic when at least one 
of these quantized momenta coincide with at least one of 
the Dirac points, defined as the zeros of the energy spec- 
trum. In the weak coupling limit U, V, V±_ <C i, t± the 
hopping Hamiltonian is sufficient for describing the band 
structure. In momentum space the energy spectrum ob- 
tained from the hopping Hamiltonian is -E'(k) = =F|ft,(k)| 
where 

h(k) = 2t cos(k x a/2)e lk y a/2V ^ + tie^ a/ ^. (5) 

Therefore the Dirac points for t = tj_ are given by 
k = (±47r/3a,0), (±27r/3a, ±27r/V3a). The Dirac points 
which corresponds to one of the quantized momenta are 
k x = ±27r/3a. This choice of quantized momenta re- 
stricts N x to being a multiple of three, i.e., a zigzag CNT 
is only metallic when N x is a multiple of three. 

The Hubbard model of a metallic CNT which is ei- 
ther armchair or zigzag may be mapped onto a two-leg 
ladder after taking a Fourier transform of the transverse 
axisi 31 ' 32 ! 33 If the annihilation operator Ci(r) is partially 
Fourier transformed in the x direction and k x — ±27r/3a, 

(,-./,) -V > :d ql (y)e^ 3 ^. (6) 



g=± 



On substitution of this Fourier transform into the hop- 



ping Hamiltonian, 

H Q =Yytd\ la {y)d q2a {y + b-) 

yqa 

-t±dl la (y)d q2a (y-b + ) + h.c] (7) 

where b± = b ± 5, b = a\/3/4, S — a/4\/3 and the Fermi 
point is kp = k y = n/2b. Similarly, the two interaction 
Hamiltonians are 

Hu =i^5ZKiT(y) n <?ii(y) + n q^(y) n ^i(y) 

yqi 

+ 4t (v)4n ^ ( 8 ) 

with n qia = d\ ia d qia and q = -q, and 
Hv = % ^2 i n qia(y)n q '2p(y + b-) 

yqq'afi 

+ 5 q , q COS x4la(?/) rf 9'l«(y)4'2^(y + b -) d q 2f}{y + &-)] 

+ w: ^2 [ n i^(y)n q '20(y - b+) 

affyqq' 

+ 6 q'qdli a (y)d q <i a (y)dl, 2fi (y - b+)d q2 p{y - &+)]. (9) 

Though the Hamiltonian now resembles a two-leg ladder 
Hamiltonian, it is not in the standard form. The two- 
leg ladder Hamiltonian usually describes hopping and 
nearest neighbour interactions both along legs and across 
rungs, and this is the type of two- leg ladder the armchair 
CNT maps onto. However, the hopping part of the our 
Hamiltonian describes a two-leg ladder with no hopping 
across the rungs and a hopping strength alternating be- 
tween t and t± along the legs, as shown in Fig. [3] Note 
that although the two legs appear to be like two inde- 
pendent chains this is not the case as they influence each 
other through slightly complicated Coulomb interactions. 

When considering weak interactions a two-leg ladder 
hopping Hamiltonian must usually be diagonalized so 
that it can written in terms of two decoupled bands. In 
our case the two legs are already decoupled in the hop- 
ping Hamiltonian making any diagonalization unneces- 
sary. Instead we can immediately make a first approxi- 
mation by linearizing the lattice fermion operators oper- 
ators about the Fermi points. This is done by expanding 
in terms of chiral fields^ 

d q i(y)/Vb ^ Rq (y)e ik ™ + ^ L q( y y~ tkFy 
d q2 (y ± b - S)/Vb ~i> Rq {y ± b)e lkF ^ 

+ xP Lq (y±b)e- ik ^ y±b \ (10) 

On substituting the chiral fields into the Hamiltonian 
we discard the rapidly varying terms, keeping in mind 
that the spatial coordinate may now be written as y — 
2mb for integral m. We take the continuum limit, which 
involves expanding ipp q (y± b) in a Taylor series about y, 
and retain the lowest order non-zero term. Finally, the 
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or '2' 



t± t± 
t t 



y + 2b 
y + b-S 
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FIG. 3: A two- leg ladder which is equivalent to a metallic 
zigzag CNT with hopping t,t±. We do not show the interac- 
tions here as they are a little complicated. The three values 
on the bottom right give general y coordinates of the lattice 
points. 



Hamiltonian density Ti may be obtained by converting 
the discrete variable y into a continuous variable so that 
Ha + Hu + Hy = fdy[Ho + Hu+Hv}/2b. The resulting 
hopping Hamiltonian density is 



(11) 



where v — tb is the Fermi velocity. 

The interaction Hamiltonians may be written in terms 
of the currents 



J qq> = Wpau^Pl'^ Jqq> = §^ko°a/3^ 



Iqq' — \4>Pqa£aplljPq' 0, Igg' — \^Pqa {^c)al3' i PPq' 



>Pq'fi 



(12) 



After combining both the on-site and nearest neighbour 
interaction Hamiltonian densities^ 



Hi =Hu + Tlv 



>er -wR -rL 

qq' qq' qq' 
R tL 



=b p , J R , J L , - 
qq' qq qq 

i fP jP jL _ to tK ~ih 

~ J qq 1 qq V 'q' J qq qq q'q' 
i l.p r-Rf jL 

'"qq'^qq' -~qq' 



, P JR]JL 
2 L 'qq' qq' qq' 



h.cl 



(13) 



where the Hermitian conjugate (h.c.) only refers to the 
final two terms which are the Umklapp interactions. In 
deriving this Hamiltonian we took the zeroth order con- 
tinuum limit of the nearest neighbour interactions. 

The coefficients f qq > and b qq i describe forward and 
backward scattering, respectively, where q = ± = 1,2. 
To avoid double counting, f qq — 0. We may also take 
u a qq = as l qq = 0. Taking Hermiticity and parity sym- 
metry into account gives 612 = 621 and /12 = /21 re- 
spectively. Therefore we have nine independent coupling 



strengths which are 



Ki = 




■3V + 3V X ), 


Ki = 




V - 


■v±), 


&?2 = 


wSu + 




6f 2 = 


m- 


2V 


-v±), 


f?2 = 




-6V + W ± ), 


/12 = 




2V 


-v±) 


u u = 


B u ~ 


2V-V X ), 










u 12 = 




-V-2V X ), 


«12 = 






(14) 



When written in the form given in Eqs. ()1 1I13P the arm- 
chair and zigzag CNT look identical, but it is important 
to note that their coupling strengths are quite different 
and so one would not expect the two nanotubes to have 
similar solutions. 

For further analysis it is convenient to bosonize the 
Hamiltonian. 29 In terms of boson fields the fermion op- 
erator is 



where P 
are 



R/L 



ipp qa = n qa e^ p ^ (15) 
±. The boson anticommutation rules 



[<pp qa (y), <t>Pq'p(y')} = iPir5 qq '5 a p sgn(y - y') 

[0Rqa(y),<l>Lq>f3(y')} =i^5 qq >8 a (16) 

and the Klein factors satisfy {n qa , Kq'p} = 28 qq '5 a p. 
From this we can define a displacement field 9 qa — 
4>Rqa - 4>Lqa and a phase field p qa = 4> Rqa + (f>Lqa- Then 
we can define a charge mode 9 qp = (9 q ^ + q i)/ \f2 and 
a spin mode 9 qrT = (9 q ^ — 9 q \)j\/2 and likewise for (p. 
Finally we define V ± = {9\ u ± 9 2v )/y/2 where v = p, a 
and similarly for ip. The Hamiltonian densities in term 
of the boson fields are 



Ul =^E^±PA±) 2 -(^±) 2 ] 



8tt *-f L 



2Tb1 2 cos </? p _ cos 6 a+ 

2 cos 9 a+ (6^1 cos 9 a - + r/^ cos p a - ) 

cos ip p - (Tb^ 2 cos 9 a _ + b^ 2 cos p a - ) 

2Yu p lx cos 9 p+ cos ip p - — 2u°2 cos 9 P+ cos 9 a+ 

cos 9 p +(uf 2 cos 9 a - + Tui 2 cos fa-)] (17) 



where A p± = b p n ±f p 2 , A a± = -(^±/f 2 ), bf 2 = b° 12 ±b p 12 
and uf 2 — u\ 2 ± u p 2 and T = ki|Ki|K2T k 21- As T 2 = 1 
we have T = ±1, although we shall set T = 1. 



III. PHASE ANALYSIS 

A. SO(8) Gross-Neveu model 

We use a well known RG treatment in order to deter- 
mine the ground state phases of our two-leg ladderj 39 i 40 i 41 
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FIG. 4: Numerical solution of the RG flow equations for the 
D-Mott phase. The cutoff is Id = 39.0/(7 and the (-axis is 
over the range A( = 0.2/17. 



FIG. 5: Numerical solution of the RG flow equations at the 
transition between the D-Mott and the S-Mott phases. The 
cutoff is Id = 35.2/(7 and the (-axis is over the range A( = 
4/(7. 



The renormalization group flow equations are equivalent 
to those obtained in Ref. once the slightly different 
definitions of Jp qq > and Ip qq >, as given in Eq. (fT2)l . are 
taken into account. To solve these flow equations we in- 
sert the ansatz 



qq' 



(18) 



where ( is the flow parameter and g" , is one of the nine 
independent coupling strengths with the initial values 
g qq i(0) given in Eq. (fl"4|) . The constants Id and G qq , are 
obtained on substituting the ansatz into the flow equa- 
tions. We solve the flow equations numerically for var- 
ious values of U, V and V± while always maintaining 

\u\>\v\,\v ± \. 

Two examples of our numerical solutions are shown 
in Figs. H] and [5] in the vicinity of I — Id- What we 
are really plotting in these nine graphs is 3(ld — 0flgg'(0- 
Fig. [4] shows that and F° 2 both flow to zero and are 
therefore negligible. Meanwhile, all other coupling con- 
stants flow to the same absolute value. As we shall show 
later, this solution of the RG flow equations describes 
the D-Mott phase. Fig. [5] describes the phase transi- 
tion between the D-Mott phase and the S-Mott phase. 
In each phase or at each phase transition the nine cou- 
pling constants flow to a specific set of values. In total 
we observe seven distinct phases, as shown in Fig. [6] for 
positive U and Fig. [7] for negative U . There are only 
four physically possible phases, the D-Mott, the S-Mott, 
a CDW and a PDW, i.e, phases obtained for positive (re- 
pulsive) interactions, though the PDW exists in a very 
narrow parameter range. The unphysical phases are the 
CCP, the D'-Mott and the S'-Mott. The defining cou- 
pling constants G qq , of each phase are given in the last 
column of Table |TJ An explanation for the names of the 
phases will be given in Sec. IIII Bl 

Most of the zigzag CNT phases (including the unphysi- 
cal ones) are also known to be armchair CNT phases. The 
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FIG. 6: Phases for |V|, \V±\<U,U> 0. 



CCP is one phase which has not been explicitly noted 
in the armchair CNT, but the SF phase in a standard 
two-leg ladde r 34 ' 35 should map to a CCP in the armchair 
CNT. One phase which can in principle exist in a zigzag 
CNT but which we do not find in our phase diagrams, 
Figs. [5] and is the FDW. Although, we do not ob- 
serve this phase for the half-filling case considered here, 
we will later show that it can be found in phase diagrams 
of doped zigzag CNT.^I 

If the coupling strengths for the D-Mott phase are sub- 
stituted into Eq. (jTTJ) it can be shown that, if we define 

(0,<p)i = (0,(f) p+ , {Q,<ph= (0,(f)a+ 

(9, cp) 3 = {9, (9, <^) 4 = 0) p _ (19) 



then 

H °=ir £KW 2 + (^a) 



n r = - 



9 
2tt 2 



^ d y (j)Rady(j)La - 4:irg ^ cos #q cos ' 



(20) 
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TABLE I: Each phase's ground state has four pinned fields and four which fluctuate rapidly. For all phases (0p_), (<Pp+) and 
(</V+) fluctuate rapidly. The fourth rapidly varying field is indicated in the table by The phase is determined by the 
coupling constants. Those coupling constants not mentioned in this table remain small so are negligible. 



phase 


<0p+> 


(0*+) 


{Oct-} 


(Vp-) 


(fa-) 


coupling strength 






D-Mott 














- 


B{ 2 = flf a = F p 2 = —B xl = E/& 


- n p 

— u 12 


= E/f a > 


S-Mott 











7T 




-Bf a = = = -Bfi = -LTfi 


— u p 

— u 12 


= U? 2 > 


D'-Mott 


7T 













B p 12 = Bf a = Ff 2 = = = - 


TT P - 

u 12 — 


-U&>0 


S'-Mott 


7T 








7T 




-B p 12 = -B? 2 = F p 2 = — -Bfj = E/& = - 




-u? 2 >o 


CCP 
















—B p 2 = Sf 2 = F12 = —Fl2 = Ull = 


U 12 


= U? 2 > 


CDW 










TV 





B? 2 = B\ 2 = F[ 2 = ~Fi2 = = 


u 12 


= U? 2 > 


FDW 


7T 













-Bf a = = = -f£ = = 


U p - 

u 12 — 


-U? 2 >0 


PDW 


7T 







7T 





S 12 = — B12 = F p 2 = —F[ 2 = U p x = 


u 12 — 


-U? 2 >0 



3 
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FIG. 7: Phases for |V|, |Vj_| < \U\, U < 0. 



where we have defined a chiral field (j>p a = (ip a + P9 a )/2 
and 4g — \G qq ,\ for G v qql ^ 0. In the semiclassical limit 
the ground state can be determined by minimizing the 
above Hamiltonian. It is not hard to see that the Hamil- 
tonian will be minimized when we simultaneously pin 
either 9 a = 2n a ir or 9 a — (2n a + l)ir for all a and integral 
n a . These pinned fields describe the gapped (or massive) 
excitations of the system as any change in a must be 
finite. The system is unaffected by any change in the 
unpinned fields, even if the changes are infinitesimally 
small, so these fields define the gapless (or massless) ex- 
citations. The ip a field is dual to 9 a , which means if one 
field is pinned the other must vary rapidly in the semi- 
classical limit i.e., the Heisenberg uncertainty principle 
comes into effect and causes one field to vary rapidly if 
its dual field is well known. 

We now show how the bosonized Hamiltonian may 
be mapped onto an SO(8) GN model. 29 Firstly we 
refermionize using 

0Pa = K Q e 1 ^, a =1,2, 3 

^ P4 = P K4 e^ 4 (21) 

where the Klein factors are k% — K21, re 2 = Kit, re 3 — i^ii, 



and K4 = K2i- We then map to Majorana fields, 

lJ->Pa = (VP2a + il]P2a-l) / V~2 (22) 

to obtain the SO (8) GN model, 

U = IvRAidyVRA - ^VLAidyVLA + gG% B G£ B (23) 

with A, B = 1, 2 . . . , 8 and currents Gp B — ijpa'Hpb for 
A ^ B, thus showing that the D-Mott phase of the two- 
leg ladder has its symmetry enlarged to SO (8). 

The Hamiltonians of the seven remaining phases in 
terms of a and ip a or rjp a may be obtained using appro- 
priate mappings from the D-Mott phase, thereby showing 
that all phases have an SO (8) symmetry, although they 
do not share the same SO (8) symmetry. If 9 a and ip a are 
as defined for the D-Mott phase in Eq. (fT!)|) then for the 
D'-Mott phase, 



a = 1 

a = 2,3,4 



(24) 



and <p® = <p a . In terms of the GN Majorana fields, 



rP' 

Vpa 




A = 1,2 
A = 3,.. 



(25) 



Therefore the two D-Mott phases share an SO (6) subal- 
gebra. Similarly for the S-Mott, 



a = 1,2,3 
a = 4 



and (fa = <p a . In terms of the GN Majorana fields, 



"Q S pa 



TIP A, 

PVpa, 



A 
A 



1,2,3,4,5,6 
7,8. 



(26) 



(27) 



Therefore the D-Mott and the S-Mott also share an SO (6) 
subalgebra, although it is different from the subalgebra 
shared by the D-Mott and the D'-Mott. All the symme- 
tries between all the phases are shown in Table HT1 
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TABLE II: Symmetries shared by different phases. The upper right corner is for the undoped case and the lower left corner is 
for the doped case. 



phase 


D-Mott/SC 


D'-Mott 


S-Mott/SC 


S'-Mott 


CCP 


FDW 


CDW 


PDW 


D-Mott/SC 




SO(6) 


SO(6) 


SO(4) 


SO(7) 


SO(5) 


SO(5) 


SO(3) 


D'-Mott 






SO(4) 


SO(6) 


SO(5) 


SO(7) 


SO(3) 


SO(5) 


S-Mott/SC 


SO(4) 






SO(6) 


SO(5) 


SO(3) 


SO(7) 


SO(5) 


S'-Mott 










SO(3) 


SO(5) 


SO(5) 


SO(7) 


CCP 


SO(5) 




SO(3) 






SO(6) 


SO(6) 


SO(4) 


FDW 














SO(4) 


SO(6) 


CDW 


SO(3) 




SO(5) 




SO(4) 






SO(6) 



B. Classification of phases 

In Sec. [J we gave a qualitative description of the eight 
phases, while in the previous section the phases were de- 
fined in terms of their pinned fields. In this section we 
relate the pinned fields to the qualitative description. We 
do this by calculating various order parameters such as 
the current, kinetic energy, electron density, supercon- 
ducting order parameter and the pair field operator. 

1. Electron density 

Each lattice site contributes one electron so the average 
electron density per lattice site is one. Deviations from 
this average value can be define by 

n(R') = ^ C L(R')c M (R') (28) 

a 

on a generic site of the carbon lattice with i — 1,2. Af- 
ter using the same mappings which provided us with the 
two-leg ladder Hamiltonian and then bosonizing the re- 
sulting equation we can represent the deviation from av- 
erage density by 

n(m) =(-l) m X]^« ? a[e- i ^ + ^ 

= -i8(-ir« 1TK2T 

x (sin \y>p- cos \y> a - cos \9 P+ cos \9 a+ 

— cos \<£p- sin \ip a - sin \0 p+ sin \0 a+ ). (29) 

The Klein factors have been simplified by using 
(k2T k it)( k it k 2|) = 1 an d r = 1. To evaluate the elec- 
tron density we consult Table H] and substitute in the 
ground state values. For all Mott phases ip p - is rapidly 
varying so n(m) = and there is no deviation from the 
average electron density of one electron per site. The 
deviation from the average electron density is also zero 
for all spin wave phases except the CDW. In the CDW 



n(m) oc (—1)™. Therefore, as shown in Fig. QJe) the 
CDW consists of pairs of electrons on alternate sites. 
Note that there are two possible forms of the CDW, one 
where all electrons are positioned on the sublattice rep- 
resented by green triangles, and the other where all elec- 
trons are positioned on the sublattice represented by red 
squares. 



2. Current between nearest neighbours 

The current between nearest neighbours on the carbon 
lattice may be defined as 

i ± (R) = -z5>L(R) C2Q (R + d)-M 

a 

ji(R) = i^[c t lQ (R)c 2Q (R + d+) -h.c] 

a 

h (R) = i [ c la (R)P2a (B. + d_) - ft.c] • (3°) 

a 

We have defined the current to always have a negative 
y component so j'j_(R) is travelling away from R along 
the perpendicular while ji.2(R) are travelling towards 
R along the zigzags. Note that when calculating order 
parameters which act between different sites we use R 
which describes all lattice sites in one sublattice of the 
CNT. In contrast, when calculating the order parameters 
which act on one site, such as the electron density, we use 
R' which describes all lattice sites in both sublattices. 
After the usual mappings we find j± = and, since R y = 
2mb, 

ji(2m) = - j 2 (2m) = «V3^ sgn(q)K qa K qa 

X ^e^^Iiia+i^Lija _|_ e i<t>Bqa -i<t>L$a j 
= «8v / 3Kl|K2t 

x (cos \y> p - cos \tp a - COS j6 p+ COS \6 a+ 

+ sin \tp P - sin \<Po- sin \9 P + sin ^6 a +)- (31) 
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This current is a function of 2m rather than m because 
we have defined the current in terms of one sublattice 
rather than the entire lattice. 

Table U shows that for all Mott phases ip p - is rapidly 
varying so ji j2 = 0. In all other phases, except the 
CCP we also 'find j 1>2 = 0. In the CCP ji = -j 2 = 
iSvjKii-KaT- This describes currents running along the 
zigzags of the CNT. Note that although the carbon lat- 
tice has a non-zero flow of current along the zigzags, all 
currents in the equivalent two-leg ladder vanish. This 
is because the total current in the y direction is zero in 
the CNT and the legs of the two-leg ladder are in the y 
direction. We have shown that there are two possibili- 
ties for the CDW (positive electron densities on one of 
two sublattices) and similarly, there are two possibilities 
for the CCP with currents flowing in one of two direc- 
tions about the CNT. Note that currents along different 
zigzags is always in the same direction, i.e., either clock- 
wise or anti-clockwise and not both in the same CNT. 

If we constructed currents similar to Eq. (|30[) for an 
armchair CNT we would find a different type of CCP 
where currents still flow around the tube, but along the 
'armchair' bonds. In the equivalent two-leg ladder we 
would find a SF phase. The equivalent two-leg ladder 
phase of the zigzag CNT's CCP cannot strictly speaking 
be classified as a SF phase as there are no true plaque- 
ttes around which current can flow. However, we could 
possibly define a plaquette in this two-leg ladder as be- 
ing between two nearest neighbours with an infinitesi- 
mally narrow width in the x direction. Then current 
flow around a plaquette is equivalent to equal currents 
flowing back and forth between two sites and cancelling 
each other out, as shown in Fig. QJg). 



3. Kinetic energy between nearest neighbours 
The kinetic energy is defined similarly to the current 
B±(R) = iJ2[c[ a (K)c 2a (K + d)+h.c] 

a 

Bi(R) = z^[c t lQ (R)c 2Q (R + d + ) + / l .c] 

a 

B 2 (R) = z^[cL(R)c 2Q (R + d_) + M. (32) 

a 

After the appropriate mappings we find Bi(2m) — 
B 2 (2m) = -B±(2m)/2 and 



') 



q.a 



= - 8«11-K2t 

x (sin \0 p+ cos 2#ct+ sin \<p p - cos \<p a - 

+ cos \0 p+ sin \0 a+ cos \<p p - sin \<£a-) (33) 

which vanishes in all phases except the PDW. In the 
PDW Si (2m) = -8ki T k 2T - Like the CDW and the CCP 



the PDW also has two possibilities. One possibility has 
positive kinetic energy on the diagonal bonds but nega- 
tive on the vertical bonds, and the other possibility has 
the signs exchanged. 



4- Current between next-nearest neighbours 

We can define six currents between next-nearest neigh- 
bours, 



iu(R) 



[c T lQ (R + a_)ci Q (R) - h.c] 



fa (R) =i t ct i« ( R ) c i« (R + a + ) — h.c] 

a 

iix(R) =*XM«( R + a-)ci«(R + a+) - h.c] 

a 

jai(R) =i^[cl a (R + d_)c 2ct (R + d) - h.c] 

a 

j 22 (R) =i£>L(R- + d)c 2Q (R + d+) - h.c] 

a 

j 2x (R) =i E [4 Q (R + d_)c2a(R + d+)-/i.c]. (34) 

a 

The first three currents act between sites of the first sub- 
lattice and the other currents act between sites of the 
second sublattice. Currents j% x and j 2x are in the x di- 
rection. All currents are defined to have negative x com- 
ponents. On expanding in terms of tpR/L we can see that 
jn = hi = -Jix = 321 = j 22 = -32x- Bosonizing gives 

jn(2m) ^v^V] sgn(q)n qa Kq a 



x e 



qa 

i<l>Rqc,+i<j>Lqa 



i<t>Rqc,-i<l>Lqa 



x (cos \<*p P - cos sin \0 p+ cos \Q a + 

— sin \tp P - sin \<p a - cos \0 p+ sin ^0 a+ ). 

(35) 

These currents vanish in all but the FDW. In the FDW 
jn(2m) = —i8y3KiiK 2 i which is equal in magnitude to 
the nearest neighbour current j\ in the CCP. Like the 
CCP all currents in the two-leg ladder vanish. 

Like the other charge density waves the FDW has two 
possibilities. The current about one sublattice can ei- 
ther be clockwise or anti-clockwise and the current about 
the other sublattice must then be either anticlockwise or 
clockwise, respectively. One could define this phase to be 
a type of CNT SF phase where the plaquettes are the tri- 
angles of either sublattice in the CNT. Both sublattices 
carry their own staggered current fluxi2& 



5. Kinetic energy between next-nearest neighbours 

The kinetic energy between next-nearest neighbours is 
similar to the current between next-nearest neighbours. 
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We define 








Bn(R) 


=*E[ c U R 

a 


+ a_)c Xa (R) 


+ h.c] 




=*Et c U R ) 


ci Q (R + a + ) 


+ h.c] 


Six (Ft) 


=*Et c U R 


+ a_)ci Q (R- 


f a + ) + h.c] 


S 2 i(R) 


=*Ei4( R 

a 


+ d_)c 2a (R- 


fd) + h.c] 


B 22 (R) 


=*E[4a(R 

a 


+ d)c 2a (R + 


d+) + h.c] 


£?2x(R) 


=^Ei4(R 


+ d_)c 2a (R- 


fd + ) + /i.c]. (36) 



Oi 



We can show that all the kinetic energies are equal, Bn — 
B\2 = Bi x = B21 = B22 = B 2x and 

B u {2m) = - i^[ K , a %e- lfe - +! *"° 

qct 

rh qa rx 'qa^ J 

= - 8ki T k 2 t 

x (sin \f p - cos g <P<r- cos cos \9 G+ 

— cos 5<y2 P - sin \f a ~ sin sin ^9 a+ ) 

(37) 

which vanishes for all phases except the CDW where 
Bn(2m) = — 8ki|«2|, which is equal to the kinetic en- 
ergy B\ 2 in the PDW. We have not included this kinetic 
energy in the CDW illustration given in Fig. [IJe) to avoid 
cluttering the picture. This kinetic energy in the CDW 
joins sites with equal electron density. 

6. Superconducting order parameter and pair field operator 

The S-SC order parameter describes pairings on the 
same site of the CNT lattice, 

As(R') = c iT (R')c ii (R / ) = E[ A «? + A i?l ( 38 ) 

where Ap 9 = %pp q ^%pp q i i is the pair field operator in leg 
q of the two-leg ladder. The D-SC order parameter de- 
scribes pairings across neighbouring sites and can also be 
written in terms of the pair field operator, 

A D± (R) =c 1T (R)c 2i (R + d) = i^[A s ,- A Lq ] 

q 

A m (R) =c 1T (R)c 2i (R + d+) 

= -*E[ A ^- A ^ W3 

A D2 (R) =c 1T (R)c 2i (R + d_) 

= -zEl A ^- A ide lW3 - (39) 
q 



Recall that R describes only one sublattice, while R' 
describes the entire CNT lattice. 

The singlet pairing between electrons in the two-leg 
ladder is defined by the sign of (ApiAp 2 ). A negative 
value indicates a (i-wave symmetry while a positive value 
indicates a s-wave symmetry. We define this symmetry 
in terms of the equivalent two-leg ladder, in analogy with 
the standard two-leg ladder^ For all four Mott phases 
the fields 9 a ± are pinned which implies that the two spin 
fields 9 qa are also pinned. Having 9 qa pinned indicates 
that excitations in the gth leg of the two-leg ladder with 
non-zero spin require energy, implying singlet pairings in 
the Mott phases. For the D-Mott and the D'-Mott phases 
(ApiAp,) < so they have d-wave symmetry. For the 
two S-Mott phases however (ApiA P2 ) > so they have 
s-wave symmetry. This explains the prefixes of the four 
Mott phases. In the four density waves (p a - is pinned 
rather than 9 a ~ so we can make no conclusion about the 
total spin in each leg. Because 9 a - is dual to ip a - it 
varies rapidly in the density waves and (ApiA^ 2 > = 0. 

From Table |T] it can be seen that the only difference be- 
tween the D-Mott (S-Mott) and the D'-Mott (S'-Mott) is 
the change in 9 P + from to ir. The field 9 p+ describes the 
charge gap and indicates that the singlet's center of mass 
is shifted by ir, or half a unit cell, when comparing the D- 
Mott (S-Mott) to the D'-Mott (S'-Mott). In the S-Mott 
the electron pairs oscillate between nearest neighbours. 
An S'-Mott is obtained from the S-Mott by shifting the 
centre of mass of each electron pair by half a unit cell 
so the S'-Mott phase must describe pairings which oscil- 
late between next-nearest neighbours. The relationship 
between the D-Mott and the D'-Mott is similar. 

We now return to the superconducting order parame- 
ters in Eqs. (|38j) and (f39|) . On bosonizing it can be seen 
that the S-SC order parameter is 

A s (m) =E«gT K 3l[ e ^ R9T+ ^ £?J + e^^+^i] 
q 

=4Ki T K U e i¥, "+ /2 

x (— cos \<f>p- sin \9 a+ sin \9 a - 

+ i sin \tp p - cos ^9 a+ cos ^9 a _)]. (40) 

The D-SC order parameters are 
q 

= - 4Ki T K U e^+ /2 

x (cos \tp p - cos \9 a + sin \9 a - 

+ isin \<Pp- sin \9 a+ cos \9 a -) (41) 
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and if j = 1,2 

1 

- ^^e^^wje'- 1 )^ 3 
=2K 1T K U e^"+ /2 

x (cos ^Vp- cos \0 a+ sin 3^0-- 
+ i sin sm 5^<r+ cos \®o- 

+ ( — 1) J V3(C0S 5<^p- COS 56*0-+ COS 5^0-- 

— ism jffto- sin 5^+ sin \6 a -)). (42) 

Because of the presence of the rapidly varying (p p + all su- 
perconducting order parameters will vanish in all phases. 



C. Phase transitions 

When numerically solving the RG flow equations we 
find seven phases defined by the coupling strengths G v , . 
The transitions between these phases are also defined 
by a unique set of coupling strengths as given in Ta- 
ble IIII1 From this table it can be seen that although 
we observed eight different transitions, there were only 
three types of transitions, defined as Gaussian, Ising and 
SO(5) GNxSO(3) WZW. The Gaussian and Ising tran- 
sitions were discussed in Ref. [29|. 

We can make several generalizations about what type 
of phase transition should exist between two particular 
phases. If two phases have the same pinned fields in the 
ground state and all but one of these fields are pinned 
to the same value then the transition between these two 
phases is Gaussian. For example, in the D-Mott phase 
6 p+ = and in the D'-Mott phase p + = ir, but all 
other pinned fields have the same values. However, if the 
coupling constants of the D-Mott <-> D'-Mott transition 
shown in Table IIIII are substituted into Eq. (JT7J) no sinu- 
soidal term involving either 9 p+ or its dual field Lp p + ap- 
pears and therefore these fields are massless at the tran- 
sition. Those fields which are pinned to the same value 
in both the D-Mott and the D'-Mott phases remain in 
the Hamiltonian so are massive. Refermionizing this new 
Hamiltonian and mapping to Majorana fermions gives an 
SO (6) GN model. In other words, this phase transition 
is described by a single gapless bosonic mode, which im- 
plies a central charge c = 1. Such a phase transition is 
Gaussian. By comparing changes in pinned fields is not 
hard to see why the first six transitions in Table HTT1 must 
all be of the same type. 

When two dual fields exchange their pinned value we 
have an Ising transition. For example, the D-Mott phase 
and the CCP have the same pinned values except in the 
D-Mott phase 6 a - = and ip a _ is rapidly varying and 
in the CCP (p a - = and # CT _ is rapidly varying. If we 
substitute the coupling constants from Table HTT1 into Eq. 
(|17p we obtain a rather complicated Hamiltonian. Wc 



can integrate out all modes which are massive at the crit- 
ical point (i.e, those terms only containing fields which 
do not change their values across the transition) to ob- 
tain an effective interaction Hamiltonian^ This effective 
Hamiltonian is refcrmionized and mapped to Majorana 
fermions and we obtain a GN model with one gapped 
Majorana fermion and one gapless. The discarded mas- 
sive modes provide six more gapped modes. Therefore, 
these types of transitions have a single gapless fermion 
indicating central charge of c = 1/2, which defines an 
Ising transition. 

We have one remaining phase transition, between the 
D-Mott phase and the PDW. This case initially appears 
to be more complicated than the Ising or Gaussian tran- 
sitions as we now have four fields which change their 
pinned value during the transition, 9 p+1 tp P -, a - and 
(fa- . Again we substitute the relevant coupling constants 
from Table HTH into Eq. (fTT)) . As several coupling con- 
stants are zero we obtain a relatively simple Hamiltonian 
and there is no real need to integrate out the massive 
modes, as was done for the Ising transition. One can 
proceed as one did for the Gaussian transition and imme- 
diately refermionize and then map to Majorana fermions. 
The resulting interaction Hamiltonian is of the form 

^/D-Mott— PDW ~ £ Or Gr (43) 

where A,B = 1,2,6,7,8, describing a SO(5) GN model. 
The remaining three Majorana fermions 773,4,5 do not ap- 
pear in the Hamiltonian so comprise a gapless SO (3) 
Wess-Zumino-Witten (WZW) models Therefore the 
phase transition between the D-Mott and the PDW is 
SO(5) GNxSO(3) WZW with a central charge c = 3/2. 
This may also be called a SU(2) 2 criticality or a C0S3/2 
phased The notation CmSn means a phase with m gap- 
less boson charge fields and n gapless boson spin fields 
while the subscript of SU(2)fc is the fc-level and is ob- 
tained from c = 3k/(2 + k). 

IV. DOPED CARBON NANOTUBE 

If the CNT is doped we move away from half-filling 
and Uf deviates from it /2b so we can no longer include 
Umklapp interactions in our Hamiltonian^ We assume 
infinitesimal doping so that the Hamiltonian densities in 
Eqs. (fTTj) and (|13p are still valid if all three Umklapp 
coupling strengths are set to zero. A RG analysis reveals 
four phases for the doped CNT as shown in Figs. [8] and 
[5J with coupling strengths given in the right column of 
Table ITVl Substituting the coupling constants into Eqs. 
(fTTj) and (fT3"j) and bosonizing the Hamiltonian gives, for 
example for the D-SC phase, 

^0 =^Yy(d y e a f + {d yV>a f) 

a 

Hi = - ^ 9y(j)Rady(f)La ~ 4# cos 6 a cos 6 b (44) 
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TABLE III: The phase transitions described here refer only to the transitions displayed in Figs. |6]and[7]and are not exhaustive. 
For example, one would expect a Gaussian transition between a CCP and a FDW and an Ising transition between a PDW and 
an S'-Mott phase, but as we did not find these transitions in our phase diagram we will not discuss them further. 

phase transition coupling strength classification 



D-Mott <- 
D-Mott <- 
S-Mott <- 
D'-Mott «- 
PDW <- 
D-Mott 
S-Mott <• 
D-Mott - 



• S-Mott 
D'-Mott 
S'-Mott 

■ S'-Mott 
CDW 
-> CCP 

♦ CDW 
-> PDW 



r 12 



F12 = B p 1± = 1/4, C/f 2 = E/f a = -Bfi = 1/2 
F? 2 = Bfi = 1/4, B p 12 = S% = -B& = 1/2 
= Bfi = 1/4, B£, = Bfa = Bfi = -1/2 
*Ta = B p n = 1/4, C/f 2 = Uf 2 = = -1/2 

= -Bfi = 1/4, Bf 2 = -*S = = 1/2 

= Bfa = £/& = [/fj = 2/5, Ff 2 = = -1/5 

Bia = f/f 2 = -t/fi = 2/5, Ff 2 = B& = -1/5 



Ff 2 = Bfa = Uf 2 = = 2/3 



Gaussian 
Gaussian 
Gaussian 
Gaussian 
Gaussian 

Ising 

Ising 

SO(5) GNxSO(3) WZW 




example, the CCP compared to the D-SC is 



for a, 6 = 2,3,4 with 9 a and ip a defined as in Eq. (|19p. 
This new Hamiltonian is similar to Eq. I|20p for the D- 
Mott phase, except that 9\ — 9 p+ and ipi = ip p+ do 
not appear. The pinned values of the D-SC phase and 
all other observed phases are easily found and given in 
Table ED 

We have shown how a D-Mott insulator may be 
mapped onto a SO (8) GN model. Similarly we can show 
that a D-SC may be mapped onto a SO (6) GN model of 
the form 



oCCP 



9 a, 


a = 


2,4 


<Pa 


a 


3 


<Pa, 


a = 


2,4 




a = 


3. 



In terms of the GN Majorana fields 



n CCP 

Vpa 




A = 6 

A = 3,4,5,7,8. 



(46) 



(47) 



H = ^RAidymA - ^VLAidyVLA + gGi B Gi B (45) 



with A, B = 3, 4 . . . , 8. In the undoped case we were able 
to find mappings from the D-Mott insulator to the seven 
other phases. Similarly the D-SC can be mapped to the 
S-SC, the CCP and the CDW in the doped case. For 



So the D-SC and the CCP share an SO (5) subalgebra. 
Symmetries between all doped phases can be found in 
Table IE 

Essentially one can think of the phases in the doped 
CNT as being formed from the combination of two of the 
undoped CNT phases. This is because in the undoped 
case each phase can be paired with another phase which 
only differs by the pinned value of 9 p+ , but in the un- 
doped case 9 p+ is not pinned. For example, the CCP of 
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TABLE IV: Pinned fields and coupling strengths of the phases in a doped carbon nanotube. 



phase 



coupling strength 



D-SC 
S-SC 
CCP 
CDW 





TV 



7T 



B p 2 — B12 — 2Ff 2 — 



D 12 



2F P 

Ar 12 



= —B12 
-B? 2 = Bf a = 2Ff 2 = 
Bf 2 = —5^2 = 2F 12 = 



= 



-25^ > 
-2Bfj > 
> 

-1B{ X > 



the doped CNT can be thought of as a combination of the 
CCP and the FDW of the undoped CNT. Similarly, the 
CDW of the doped CNT is similar to a combination of 
CDW and PDW of the undoped CNT and the D(S)-SC 
is like a combination of the D(S)-Mott and D'(S')-Mott 
phases. 

When we evaluate the order parameters of the doped 
CNT ground state phases we find that they are compara- 
ble to the order parameters of the undoped CNT phases. 
Using the formulas derived in Sec. IIIIBI we find that in 
the CCP of the doped CNT 



ji(2m) = i8v / 3Ki|K2| cos \0 p + 



jn(2m) 



k 2 t sin 2^ P + 



(48) 



with ji = j 2 and j n = ji2 = -jix = hx = hi = ~hx 
and all other order parameters vanish. Therefore the 
CCP of the doped CNT contains both nearest neighbour 
currents like the CCP of the undoped CNT and next- 
nearest neighbour currents like the FDW of the undoped 
CNT. Unlike the undoped case 9 p+ is not pinned so these 
currents are not set to some fixed value. Similarly, in the 
CDW of the doped CNT 



B~i (2m) = — 8ki|K2| sm \®p+ 
-Bn(2m) = — 8/ci|K2t cos \& P + 

B12 



(49) 



B\ x 



with Bi = B 2 = Bj_/2 and B u 
B 21 = B 22 = B 2x and n(m) oc (-l) m cos \0 p+ . All 
other order parameters are zero in the CDW. In the 
D-SC and S-SC all currents and kinetic energies vanish 
and the deviation from average electron density is zero. 
The superconducting order parameters in the D-SC are 
As (to) = Afli(2m) = and 

A m (2m) = -Au 2 (2m) = -2V3/ii T K U e' i¥, '+ /2 (50) 

which is not fixed because ip p + is not pinned. In the 
S-SC, 



As(to) = 4iKi T K U e 4V "+ /2 



(51) 



while all the D-SC order parameters are zero. 

The phase transitions in the doped nanotube are simi- 
lar to the transitions in the undoped system. A list of all 
the transitions and their critical values is given in Table 



IVl As before, those transitions which require one field's 
pinned value to change from to n are Gaussian. The 
massive modes at a Gaussian transition can be shown 
to be described by a SO (4) GN model. Those transi- 
tions in which a is replaced with its dual field ip a are 
Ising. The massive modes at the Ising transition are de- 
scribed by an SO (5) GN model. The transitions which 
are neither Gaussian nor Ising are D-SC<->CDW and S- 
SC^CCP. These two transitions can be understood by 
substituting their coupling constants into the Hamilto- 
nian density and then mapping to a GN model, as was 
done for the D-Mott<->PDW transition in the undoped 
case. We find that the Hamiltonian of both transitions 
resembles an SO (3) GN model, making the transition 
SO(3) GNxSO(3) WZW, or C0S3/2. 



SOLITON SPECTRUM AND QUANTUM 
NUMBERS 



In this section we present some general features of the 
SO (2 AO GN model with special reference to SO (6). De- 
tailed derivations and analyses can be found in Refs. 
USUI The SO(2A0 GN model with integer N is in- 
tegrable so the soliton excitation spectrum can be cal- 
culated exactly. These soliton excitations can be related 
directly to the CNT by considering quantum numbers 
such as charge and spin. The SO (8) case has been dis- 
cussed in Ref. [29l 



A. Excitation spectrum 

A semiclassical analysis of any Hamiltonian of the form 
given in Eq. (|44p determines the ground state to be where 
all 6 a = 2n a TT or all 6 a = (2n a + 1)tt. Therefore we have 
several possible solutions for the set of all 6 a for any given 
ground state. The system can move between any two so- 
lutions (in the one ground state) by emitting or absorb- 
ing a particle comprised of one or several solitons. The 
properties of these particles are defined by the changes 
in 6 a . 

For SO(2A jr ) a — 1,2, ... N we consider the change in 



14 



TABLE V: The phase transitions in the doped case 



phase transition 



coupling strength 



classification 



D-SC <-> S-SC 
CDW <-> CCP 
D-SC <-> CCP 
S-SC <-> CDW 
D-SC <-> CDW 
S-SC <-> CCP 



*ia = -#11 = = -Bfi = 1/3, Sf 2 = 2/3 
= -Bfi = *S = Bf a = 1/3, flf a = -2/3 
^i P 2 = -Sfi = 1, Bi P 2 = 2 



^12 



1 R p 

■ D 12 



Gaussian 
Gaussian 

Ising 

Ising 

SO(3) GNxSO(3) WZW 
SO(3) GNxSO(3) WZW 



9 a from y = — oo to y = oo, 

/oo 
dt/0 w o = 2nN a . 
-oo 

where we define the charges N a as 

V„ / dy^ a , a=l,2,...N. 



by e n (q) = \/m 2 n + q 2 ■ 



(52) 



(53) 



The N dimensional vector (N\, iVjj, . . . , iVjy) defines a 
soliton. The simplest solitons have only one # a chang- 
ing by 2n over the y range and all the others remain- 
ing constant. This corresponds to one N a = ±1 and 
all others zero. These are defined as the fundamental 
or elementary particles and there will be 2N of them. 
E.g., for SO(6) the fundamental particles are (±1,0,0), 
(0, ±1, 0) and (0, 0, ±1). Another type of soliton changes 
9 a = 2n a TT to 6 a = (2n a ± 1)tt for all a over the 
y range so N a = ±1/2 for all a. Such a soliton is 
called a kink. The kinks have no zero charges and each 
N a has two possible values and therefore there are 2 N 
kinks. The fundamental fermions and kinks are collec- 
tively known as solitons. Thus SO (8) has a total of 
24 solitons (16 kinks plus eight fundamental particles), 
while SO (6) has 14. The even kinks (or simply kinks) 
are defined to have an even number of positive charges, 
while the odd kinks (or antikinks) have an odd number 
of positive charges. For SO(6) the kinks are (—1, 1, l)/2, 
(1,-1, l)/2, (-1,-1, -l)/2, (l,l,-l)/2, and the an- 
tikinks are (l,l,l)/2, (-1,-1, l)/2, (1,-1, -l)/2 and 
(— 1,1,— 1)/2. Each fundamental particle may be con- 
structed from a kink-antikink pair if N is even, but if 
N is odd the fundamental particles are constructed from 
kink-kink or antikink-antikink pairs. Additional bound 
states may be constructed from other kink and funda- 
mental particle combinations. 

The masses of all particles constructed from kinks sat- 
isfy 



2msm[niT/{2(N - 1))], n = 1,2, . . . , N - 2 (54) 



where m is the mass of a kink. The fundamental 
particles are n = 1 and higher values of n describe 
other bound states. The energy dispersion is given 



e N ^(q) = 2^/m 2 + g 2 /4 



For n = N — 1 we obtain 
= e c (q) above which is a con- 
tinuum of scattering states. In the case of SO(6) the 
fundamental fermions have mass mi = \/2m which gives 
ei(q) — y // 2m 2 + q 2 and there are no higher bound states 
as 62 = e c - The highest bound states of the SO (8) case 
have rri2 = V3m and are constructed from kink-kink or 
antikink-antikink pairs. SO (8) is special as mi = m, 
hinting at the 'triality' of this group, i.e, an additional 
threefold symmetry between kinks, antikinks and funda- 
mental particles. Triality is also indicated by the kinks' 
and fundamental particles' distance from the origin, be- 
ing equal to unity for all three of these solitons in SO (8). 
For SO(2iV) the vectors (N 1: N 2 ,..., N N ) of the funda- 
mental particles are always a distance of unity from the 
origin, but for the kinks and antikinks they are a distance 
of v^/V /2. Therefore for N < 3 the kinks are closer to the 
origin than the fundamental particles, but for N > 5 they 
are further away. The decreasing distance from the ori- 
gin of the kink particles as N reduces below 4 implies an 
increasing instability in the bound states and fundamen- 
tal fermions. In fact, for N = 2 there is no fundamental 
particle spectrum as e\ = e e i 42 i 43 



B. Quantum numbers 

Various quantum numbers can be defined by the 
fermion operators ipPja arid may be re-expressed in terms 
of one of the bosonic fields 8 a or ip a ^ The electronic 
charge in the two-leg ladder or CNT is defined as 



Pja 

The spin is defined by 



Pja 

from which we obtain 

S z = A6 a+ /2ir. 



(55) 



(56) 



(57) 
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The relative z-component spin in different bands is 



VI. CONCLUSION 



s i2 = / dy^[i)\, la a^ l3 i>pi a -i) j P2a <T% fj 4>p2 a \/2 



Pa 



A0 CT _/27T, 



(58) 



while the relative z-component vector chirality in differ- 
ent bands is 

P12 = / dy^P[^ Pla a z af} ^ Pla -4 2 X^P2a]/2 



= A^ ct _/2tt, 
and the relative band chirality is 



(59) 



P12 = / dy^2,P{^ Pla i)p la - ipp 2oi ipP2 a ] 



A<p p -/2ir. 



(60) 



For all the undoped CNT Mott phases A6q = A9 p+ , 
A6 2 = A6 a+ , A6» 3 = A6» CT _ and A6> 4 = A^ p _, 
which after using Eq. |J55J gives (N t , N 2 , N 3 , i\T 4 ) = 
(Q/2,S z ,,Sf 2 ,P 12 ). As A# a = 2n Q ^,(2n a + 1)tt for in- 
tegral n a we have N a — n ai n a + |, showing that the 
kinks span all possible particles. For the undoped CNT 
density wave phases only N3 is different when compared 
to the Mott phases. In this case we have A# 3 = A<y9 CT _ 
so (JVi, JV 2j JV 3 , N 4 ) = (Q/2, S z , Pf 2 , P 12 ). 

Doping breaks the large SO (8) symmetry but despite 
this we may still represent the SO (6) solitons as we rep- 
resented the SO(8) solitons. The SO(6) GN model docs 
not contain 9 p+ so Q is not a good quantum number 
and the solitons should be defined by (S z , Sf 2 , P12) i n 
the superconducting sates or (S z , P* 2 , P\ 2 ) in the den- 
sity waves. However, this GN model does not describe 
the full Hamiltonian. To obtain the full Hamiltonian we 
must include a chemical potential term to the Hamilto- 
nian, H = Hqn — nQ, where [i is the chemical potential. 
Therefore Q is still a good quantum number of the doped 
Hamiltonian, even though the symmetry between Q and 
the SO (6) generators is broken. So, as in the undoped 
CNT phase, {N u N 2 , N 3 , N 4 ) = {Q/2,S z ,Sf 2 ,P 12 ) for 
the doped CNT Mott phases and (N u N 2 , N 3 , N 4 ) = 
{Q/2,S z ,Pl 2 ,P 12 ) for the doped CNT density wave 
phases. 

Other quantum numbers can also be derived, although 
they will not be simply related to the N a charges. For 
example, the momentum in the y direction (along both 
legs) may be written as 

P V = J d yJ2 kF i^h^Rj« -^Lja^Lja) 
jot 

= [(k F1 + k F2 )Aip p+ + (k F1 - k F2 )Atp p -}/2ir. (61) 

At half-filling the energy dispersion is degenerate and 
kpx = kp 2 = tt/2 so P y = Aip p+ /2. If the system is 
doped we move away from half- filling so kp\ 7^ k F2 . As 
(^p+ is not defined in the doped case with SO (6) symme- 
try P y is not a well defined quantum number. 



The metallic zigzag CNT with weak coupling and 
short-range interactions has a complex phase diagram. 
We classify all ground states as well as all phase tran- 
sitions for both the undoped case and the doped case 
within the parameter range |V|, ]Vj_| < \U\. We ob- 
tain these results by exploiting the relationship between 
a zigzag CNT and an unusual form of a two-leg ladder. 
Once the CNT model has been mapped to this two-leg 
ladder, well established RG and bosonizations techniques 
are applied to reveal the phase diagram. 

Previous studies on both doped and undoped armchair 
CNT and their equivalent two-leg ladders have found 
phases similar to those found here. Though we did ex- 
pect to find similar phases in these two CNT, we did not 
expect to find similar phase diagrams as the initial con- 
ditions of the RG equations for the armchair and zigzag 
CNT are very different (Eq. (fTi]) for the zigzag case). 
Surprisingly, despite very different initial conditions, the 
repulsive interactions part of our zigzag CNT phase dia- 
gram is remarkably similar to the armchair CNT phase 
diagram (with interactions U, V,V± > and t = but 
with the V and V± axes exchanged and the latter axis 
rescaled4^ We find that this is true for both the doped 
and undoped cases. In other words, the positive V and 
V± segment of Figs. |6] and [8] is similar to the phase di- 
agram of the armchair CNT with the same interactions 
and hopping but with V and V± exchanged. This is quite 
likely related to relative rotations of the graphene lattices 
in the two CNT as in the zigzag CNT V describes inter- 
actions in the longitudinal direction while V± describes 
interactions around the circumference, yet in the arm- 
chair CNT V± is around the circumference and V is in 
the longitudinal direction. 

Once the ground states have been established a 
refermionization of their associated Hamiltonians reveal 
that they can all be mapped onto a SO(2N) GN model, 
where N = 4 for the undoped case and N = 3 for the 
doped case. In both cases the symmetry is much larger 
than the symmetry of the initial Hamiltonian so both are 
examples of dynamical symmetry enlargement. Though 
all phases in the undoped case have an SO (8) symmetry, 
they do not all share the same subalgebra. Similarly, in 
the doped case no two phases share the same SO (6) sub- 
algebra. The subalgebra shared by two phases can hint at 
the symmetry of the phase transition between these two 
phases. We observed that if two SO(2N) phases share a 
SO(2N-l) subalgebra then the phase transition between 
these two phases is Ising and at the transition the massive 
fields are defined by a SO(2N-l) GN model. If instead the 
two phases share a SO(2N-2) subalgebra then the phase 
transition is Gaussian with the massive fields defined by a 
SO(2N-2) GN model. In both these cases the fields which 
are pinned to equal values on either side of the transition 
remain massive at the phase transition, while the fields 
which change their pinned values become massless at the 
phase transition. This is why the shared symmetry of the 
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two phases is equivalent to the symmetry of the phase 
transition. We cannot, however, always expect this to be 
the case. The D-Mott phase and the PDW phase share 
an SO (3) symmetry and yet at the phase transition the 
massive modes are described by a SO(5) GN model. This 
phase transition is distinctly different from the Gaussian 
and Ising transitions because in this case it is the fields 
which change their values across the transition which re- 
main massive while the fields which do not change their 
values becomes massless. 
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